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Abstract 

The summary of the results of our next-to-next-to-leading fits of the Teva- 
tron experimental data for xF% structure function of the vN deep-inelastic 
scattering is given. The special attention is paid to the extraction of twist-4 
contributions and demonstration of the interplay between these effects and 
higher order perturbative QCD corrections. The factorization and renormal- 
ization scale uncertainties of the results obtained are analysed. 
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1. The study of deep- inelastic scattering (DIS) processes has a rather long and 
inspiring history. One of the first realizations that the analysis of vN DIS could 
play an important role in investigations of the properties of the nucleon came in 
Ref. Q. The fundamental concept of scaling of DIS structure functions (SFs) H has 
lead to many subsequent investigations. Other important stages in the development 
of both theoretical and experimental studies of various characteristics of DIS pro- 
cesses in this productive period were reviewed in detail recently ||. In particular, it 
was stressed that after the experimental confirmation of scaling and indications of 
the existence of point-like constituents of the nucleon, the more rigorous theoretical 
explanation of the behaviour of DIS form factors came onto the agenda. A series of 
works by N. N. Bogolyubov and coauthors [Q, were devoted to the development of 
the new method, which made it possible to analyse the asymptotics of the form fac- 
tors of eiV DIS using the Jost-Lehmann-Dyson integral representation, and explain 
the property of scaling (or as called it by the authors of Ref. Q "automodelling" ) 
behaviour of the corresponding SFs in the framework of general principles of local 
quantum field theory 0. 

We now know that this property is true only in the asymptotic regime and 
that it is violated within the framework of QCD (see e.g. the extensive discussions 
in a number of books on the subject 0). Indeed, the theory of QCD predicts 
that scaling or automodelling behaviour of SFs is violated by the logarithmically 
decreasing perturbative QCD contributions to the leading twist operators. However, 
in the intermediate and low Q 2 regime the higher twist operators, which give rise 
to scaling violations of the form l/Q 2 , 1/Q 4 , etc., might also be important |7|, j8j. 
Indeed, the NLO DGLAP fits [| of the BCDMS data of DIS of charged leptons on 
nucleons [[HJ and reanalysed SLAC eiV data resulted in the detection of the 
signals from the twist-4 contributions. 

During the last few years there has been considerable progress in modeling these 
effects with the help of the infrared renormalon (IRR) approach (for the review see 
Ref.£E2|) and the dispersive method |13| (see also Ref-HHp. Using these methods 
the authors of Ref. |15| explained the behaviour of the twist-4 contributions to the 
F 2 SF observed in Ref. || and constructed a model for the similar power-suppressed 
corrections to xF% SF. In view of this it became important to check the predictions 
of Ref. fl5| and to study the possibility of extracting higher-twist contributions from 
the new more precise experimental data for vN DIS, obtained by the CCFR collab- 
oration at Fermilab Tevatron fl6| , and also to exploit the considerable progress in 
calculations of the perturbative QCD corrections to characteristics of DIS, achieved 
in the last decade. 

Indeed, the analytic expressions for the next-to-next-to-leading order (NNLO) 
perturbative QCD corrections to the coefficient functions of SFs F 2 and xF 3 
[ill are now known. Moreover, the expressions for the NNLO corrections to the 
anomalous dimensions of non-singlet (NS) even Mellin moments of F 2 SF with 
n = 2,4,6,8, 10 and for the N 3 LO corrections to the coefficient functions of these 
moments are also available [O]. In this report we will summarize the results of 
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the series of the works of Refs.pO||-|P4li devoted to the analysis of the CCFR data 
at NNLO, which has the aim to determine the NNLO value of the QCD coupling 
constant a s (Mz) and to extract the effects of the twist-4 contributions to SF xF 3 
21, 22 1 . In particular, we will concentrate on the discussion of the factorization and 



renormalization scale uncertainties of the results obtained. 

2. Our analysis of Refs.|20||- |22|| is based on reconstruction of the SF xF% from 



its Mellin moments M n (Q 2 ) = J x n ~ l Fz{x,Q 2 )dx using the Jacobi polynomials 
method, proposed in Ref. 
this framework one has: 



and further developed in the works of Ref. |24| . Within 
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where O"" 9 are the Jacobi polynomials, ^(a, (3) are combinatorial coefficients given 
in terms of Euler T-functions of the a and (3 weight parameters. In view of the rea- 
sons discussed in Ref.[^2|, they were fixed to 0.7 and 3 respectively. The QCD 
evolution of the moments is defined by the solution of the corresponding renormal- 
ization group equation: 
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The QCD running coupling constant enters this equation through A S {Q 2 ) = a s (Q 2 ) / (Att) 
and is defined as the expansion in terms of inverse powers of/n(Q 2 /Ag/). For 
the initial scale Qq, from which the evolution is started, the moments in Eq.(2) 
were parametrized as M n (Ql) = ^ x n ~ 2 A(Ql)x b ^o){\ — x) c ^o){\ + >y(Q 2 )x)dx. 
In the process of our analysis we took into account both target mass corrections 
and twist-4 contributions. The latter were modeled using the IRR approach as 
M* RR = C(n)M n (Q 2 )A' 2 /Q 2 |TJj and by adding into the r.h.s. of Eq.(l) the term 
h(x)/Q 2 with h(x) considered as a free parameter for each x-bin of the experimental 
data. 

For arbitrary factorization and renormalization scales the NNLO expression for 
the NS Mellin moments reads: 



M n (Q 2 ) ~ (A s (Q 2 k F )) a x AD(n,A s (Q 2 k F )) x cP s (A s (Q 2 k R )) 



(3) 



where a = 7^/(2/%), AD = 1 + 



A s (Q 2 k F ) + 



q(n) + p(n)(a + l)kf + 



(f3 1 /(3 )akf + a(a+l)(kf) 2 /2 



p(n) + ak F 

A 2 s (Q 2 k F ) and = l+C^{n)A s {Q 2 k R )+ C&(n) + 

C^\n)k R A 2 (Q 2 kn). Here 7^5, /3q and (3\ are the scheme-independent coefficients of 

the anomalous dimension function ^ns{ x ) an d QCD /3-function /3(x), p{n) and q{n)- 
terms are expressed through the NLO and NNLO coefficients of 7at,s(x) and f3{x) via 



equations, given in Refs. J2(], |22[ . Within the MS-like schemes the factorization and 
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renormalization scale ambiguities are parameterized by the terms k[ = ^oln^kp) 
and hi = f3oln(kR), where kp (&r) is the ratio of the factorization (renormalization) 



scale and the scale of the MS'-scheme. Following the analysis of Ref.|25[ we take 
kn = kp = k, fixing identically the factorization scale and the renormalization scale. 
We performed our fits for the case of k = 1 (namely, in the pure MS'-scheme) and 
then determine the scale uncertainties of A^-, the twist-4 parameter A 2 and the 
x-shape of h(x) by choosing k — 1/4 and k = 4 and repeating the fits for these two 
cases. 

3. In the process of our analysis of CCFR'97 data we applied the same kinematic 
cuts as in Ref.fl6j, namely Q 2 > 5 GeV 2 , x < 0.7 and W 2 > 10 GeV 2 . We started 
the QCD evolution from the initial scale Qq = 20 GeV 2 , which we consider as more 
appropriate from the point of view of stability of the NLO and NNLO results for 



due to variation of the initial scale pifl . In order to estimate the uncertainties of the 
NNLO results, we also performed the N 3 LO fits with the help of the expanded Pade 
approximations technique (for the detailed discussions see Ref. p2|| ). The results are 
presented in Table 1. 





A^ (MeV) 


A' 2 (GeV 2 ) 


X 2 /points 


LO 


264±37 




113.1/86 




433±53 


-0.33±0.06 


83.1/86 




331±162 


h(x) in Fig.l 


66.3/86 


NLO 


339±42 




87.6/86 




369±39 


-0.12±0.06 


82.3/86 




440±183 


h(x) in Fig.l 


65.8/86 


NNLO 


326±35 




77.0/86 




327±35 


-0.01±0.05 


76.9/86 




372±133 


h(x) in Fig.l 


65.0/86 


N 3 LO 


332±28 




76.9/86 


Pade 


333±27 


-0.04±0.05 


76.3/86 




371±127 


h(x) in Fig.l 


64.8/86 



Table 1. The results of the fits of CCFR'97 data with the cut Q 2 > 5 GeV 2 . 

At NLO the value for A^, is in good agreement with the NLO result A^, = 
337 ± 28 MeV, obtained by the CCFR collaboration with the help of DGLAP 
NLO analysis of both F2 and xF% SFs data in the case when HT-corrections were 
neglected [0. The obtained NLO value of the IRR-model parameter A 2 is in 
agreement with the estimates of Ref.flT5| and of Ref.p6| especially. However, at 



NNLO a significant decrease of the magnitude of the parameter A 2 is observed. In 
view of this the results for A^- obtained at the NNLO without HT corrections and 
with IRR-model of twist-4 term almost coincide. A similar tendency was observed in 
the process of the N 3 LO Pade fits. To study this feature in more detail we extracted 
the x-shape of the model- independent function h(x) (see Fig.l) and analysed the 



4 



factorization/renormalization scale uncertainties of the outcomes of our fits [22 
The corresponding results are presented in Table 2 where is defined as 
A^^(k) — A^L(/c = 1). The related x-shapes of h(x) are presented in Fig. 2. 

1V1 O lvl o 



Order 


k 


A k (MeV) 


A 2 (GeV 2 ) 


X 2 /points 


NLO 


4 


116 




99.1/86 




4 


213 


-0.22±0.006 


84.2/86 




1/4 


-61 




80.4/86 




1/4 


-99 


+0.02±0.005 


80.2/86 


NNLO 


4 


35 




83.5/86 




4 


66 


-0.11±0.06 


83.5/86 




1/4 


-51 




87.3/86 




1/4 


-45 


+0.09±0.05 


84.5/86 



Table 2. The results of NLO and NNLO fits of CCFR'97 data for different values 
of factorization/renormalization scales. 



4. We will concentrate first on discussing the presented behaviour of the twist-4 
parameter h(x) of xF 3 SF, presented in Figs. 1,2. In the case of k — 1, namely in 
the pure MS'-scheme, x-shape of h(x) obtained from the LO and NLO analysis of 



Refs. [pi| , |22j[ is in agreement with the IRR-model predictions of Ref.||15||. Note also 
that the combination of the quark counting rules |27j with the results of Ref.0 
predict the following x-form of h(x): h(x) ~ A 2 (l — x) 2 . Taking into account the 
negative values of A 2 , obtained in the process of our LO and NLO fits (see Table 
1), we conclude that the related behaviour of h(x) is in qualitative agreement with 
these predictions. 

At the NNLO the situation is more intriguing. Indeed, though a certain in- 
dication of the twist-4 term survives even at this level, the NNLO part of Fig.l 
demonstrates that its extracted x-shape starts to deviate both from the IRR predic- 
tion of Ref . |15| and from the quark-parton model picture, mentioned above. Notice 
also that within the statistical error bars the NNLO value of A 2 is indistinguish- 
able from zero. These conclusions are confirmed by the studies of the factoriza- 
tion/renormalization scale dependence of the NLO and NNLO outcomes of the fits 

Indeed, it is known that the variation of the related scales is simulating in part 
the effects of the higher-order perturbative QCD corrections. In view of this the 
NLO (NNLO) results, obtained in the case of k = 1/4 (see Table 2 and Fig. 2 in 
particular), are almost identical to the NNLO (Pade motivated N 3 LO) extractions 
of h(x) and of the IRR model parameter A 2 from the fits with k — 1 (see Fig.l and 
Table 1). Thus, we conclude, that as the result of analysis of the CCFR'97 data the 
NNLO and beyond we observe the minimization of the twist-4 contributions to xF 3 
SF. This feature is related to the interplay between NNLO perturbative QCD and 
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Figure 1: extracted from CCFR'97 data 
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Figure 2: Scale dependence of 
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twist-4 1/Q 2 corrections. The recent studies of the scale-dependence of the NLO 
DGLAP extraction of the twist-4 terms from different recent DIS experimental data 
f28fl are supporting the foundations of Refs.^1], p^j . This means that the higher- twist 
parameters cannot be defined independently of the effects of perturbation theory 
and that the NNLO corrections can mimick the contributions of higher twists |29| 



provided the experimental data is not precise enough for the clear separation of the 
nonperturbative from perturbative effects. Thus, it is highly desirable to have new 
experimental data for xF% SF, which are more precise than the ones given by the 
CCFR collaboration. 

In conclusion we present also the NLO and NNLO values of a s (Mz), obtained by 
us in Ref.p^j from the fits of CCFR'97 data for xF 3 SF with twist-4 terms modeled 
through the IRR approach : 

NLO a s (M z ) = 0.120 ± 0.003(stat) ±0.005(syst)±%%$ (4) 
NNLO a s (M z ) = 0.118 ± 0.003(stat) ± 0.005(syst) ± 0.003 

The systematical uncertainties in these results are determined by the systematical 
uncertainties of the CCFR'97 data and the theoretical errors are fixed from the 
numbers for (see Table 2), which reflect the factorization/renormalization scale 
uncertainties of the values of A^=-. The incorporation into the MS-matching formula 



|3H] of the proposals of Ref.[[lT]] for estimates of the ambiguities due to smooth 
transition to the world with / = 5 numbers of active flavours was also taken into 
account. The theoretical uncertainties presented are in agreement with the ones, 
obtained in Ref . |[25|| , while the NNLO value of a s (Mz) is in agreement with another 
NNLO result a s (Mz) = 0.1172 ± 0.0024, which was obtained from the analysis of 
SLAC, BCDMS, E665 and HERA data for F 2 SF with the help of the Bernstein 
polynomial technique |32j| . It might be of interest to verify the theoretical errors 



of these two available phenomenological NNLO analysis using different variants of 
fixing scheme-dependence ambiguities. The first steps towards the analysis of this 



problem are already made [33 
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